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How much multifractality is included in monofractal 

signals? 
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Abstract 

We investigate the presence of residual multifractal background for monofractal signals 
which appears due to the finite length of the signals and (or) due to the long memory 
the signals reveal. This phenomenon is investigated numerically within the multifractal 
detrended fluctuation analysis (MF-DFA) for artificially generated data. Next, the ana- 
lytical formulas enabling to describe the multifractal content in such signals are provided. 
Final results are shown in the frequently used generalized Hurst exponent h{(i) multifrac- 
tal scenario and are presented as a function of time series length L and the autocorrelation 
exponent value 7. The multifractal spectrum (a, /(a)) approach is also discussed. The 
obtained results may be significant in any practical application of multifractality, including 
financial data analysis, because the 'true' multifractal effect should be clearly separated 
from the so called 'multifractal noise'. Examples from finance in this context are given. 
The provided formulas may help to decide whether we do deal in particular case with 
the signal of real multifractal origin. They also push further findings already existing in 
literature. 
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1 Introduction 



Multifractality [HElEllllElElE] is the property of complex and composite systems that has 
been attracting more and more attention in recent years. The practical fruits of multifractality 
are not precisely known yet but at least in the case of financial markets some interesting features 
of this phenomenon were shown (see e.g. [H El lIOl IHl US [131 El [El [SI [H]) that rise hopes for 
future applications. Since the paper by Kantelhardt et. al. [IS] we know that multifractality may 
result not only from long-range correlations but also from fat tails in probability distributions 
(PDF) of investigated data. Normally, one expects multifractality in time series as a result 
of different form of autocorrelations appearing at various time scales. However, we always 
deal in practise with finite samples of data collected in time series of given length. In such a 
case, multifractality may appear even if no difference in autocorrelation properties exists for 
various time scales. It is because large fluctuations cannot be detected as frequent as small 
fluctuations in finite samples of data with long memory - mainly due to the insufficient data 
statistics. In other words, large fluctuations are not able to be formed in small samples of 
data, contrary to small fluctuations. Therefore, one gets in the case of shorter time series the 
multifractal property which itself is not programmed to be multifractal in a sense of different 
autocorrelation properties at various scales. The latter multifractality, related to variety of 
autocorrelations, is more substantial and has to be somehow separated from the former one 
which we shall call the multifractal background or residual noise further on. The preliminary 
analysis of this problem had already been made in [191 12D] • Our goal is to describe the presence 
of multifractal background quantitatively for time series with long memory induced by explicit 
form of autocorrelations in data. Our approach is directly based on Fourier filtering method 
(FFM) [2T] and differs therefore from other approaches, where the long-memory effect was 
inserted into a signal by the particular choice of power spectrum (amplitude adjusted Fourier 
transform method [22]) [19] or log- normal cascade implementation [20], instead of direct shaping 
the artificial data with autocorrelation exponent 7 discussed further on. 

We use the multifractal detrended fiuctuation analysis (MF-DFA) [4J as the commonly 
accepted technique to find multifractal properties of time series. This method is described 
elsewhere (see e.g. [H |T6l |171 [18]) so we will not recall it in details here. To keep the standard 
notation we will note the g-deformed fiuctuation of the time series signal around its local trend 
(assumed linear in our approach) in a time window of size r as Fq{T). Usually, the multifractal 
properties are presented as the multifractal spectrum (a, /(a)) [7], called sometimes also Holder 
description. Equivalently, in the Hurst language, one can consider the spread of generalized 
Hurst exponents h{q) [18], calculated within MF-DFA for g — deformed fiuctuations -Fg(r) from 
the power law: 

F,{t) ~ T^^"^ (1) 

Both descriptions are linked together via relations [231 [21] 

a = h{q) + qh'{q), f{a) = q{a - h{q)) + 1 (2) 

We start in this article with the generalized Hurst exponent h{q) description of multifractal- 
ity. Our results are then easy translated into Holder language with the use of Eq.(2). Finally, 
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we compare our findings with the properties of real data from financial market, commonly 
believed to exhibit multifractal features. 

2 Generalized Hurst exponents for finite monofractal 
signals 

Our aim is to evaluate the multifractal effect in finite artificial signals of various lengths for 
given constant level of persistency, i.e. that autocorrelations are not changing with the time 
scale. Such signals are built by us within FFM [21]. The level of autocorrelations is modulated 
by the proper choice of scaling exponent 7 responsible for the magnitude of autocorrelation 
function C(r). The latter one satisfies for stationary series with long memory the known power 
law: 

C(r) = (Ax(t)Ax(t + r)) ~ r"^ (3) 

where Ax{t) = x{t+l)—x{t) are increments of time series, r is the time-lag between observations 
and the average () is taken over all data in the series. 

The 7 scaling exponent may be linked to the Hurst exponent H [25] by the formula [26] : 

7 = 2-2/7. (4) 

In the quantitative analysis of residual multifractality left in monofractal finite signals, 
we concerned the ensembles of numerically generated time series of length L = 2", (ra = 
9, 10, ... , 20) with the pre-assumed autocorrelation exponent value 7 = 0.1, 0.2, . . . , 0.9, 1.0, 
each containing 10^ independent realizations. Thus, the spread of 7 exponents covers the range 
1/2 < H < 1. The obtained quantities have been averaged over such statistical ensemble with 
given L and 7 as input parameters. 

First, we examined the FFM procedure for time series generation, in order to check its 
accuracy towards replication of the pre-assumed autocorrelation properties coming out from 
the particular choice of 7 exponent as input. Fig.l demonstrates its efficiency. It is seen that 
the power law in Eq.(3) is reproduced very well even for large time-lags. Moreover, a coincidence 
between input and output 7's is also satisfactory. The length of generated data-samples was 
chosen as powers of 2 to improve performance of fast Fourier transform algorithm. 

The next problem we had to examine, is the performance of MF-DFA technique which 
strictly depends on the power law scaling between q — deformed fluctuations Fgi^r) and the box 
size r (see Eq.(l)). An exact extraction of the generalized Hurst exponent h{q) is then possible 
only for well determined scaling range in the fitting procedure logFq(r) vs logr. Fig.2 clearly 
shows the expected power law dependence for various lengths of the signal L and for different 
values of deformation parameter q. The latter one was uniformly distributed in the range from 
—15 to +15. These plots justify the scaling range from r = 10 till r = L/4 which was chosen 
by us to be used further on. 

To determine quantitatively the amount of multifractal residual noise present in given time 
series, the edge values of h{q) function were investigated for them. Let us introduce the new 
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parameter Ah defined as the difference between the two asymptotic hmits 

A/i = lim hiq) — lim hiq) (5) 

and assume for numerical reasons that such asymptotic hmits are reached aheady at g = ±15. 
Such assumption is justified in Fig.3, where plots for h{q) are shown for L = 2^^, 2^" and for 
7 = 0.1, 0.5, 0.9 respectively. 

Generally, we may expect that Ah is the unknown function of L and 7. The form of Ah{L, 7) 
dependence is thus a crucial problem. To simplify it, one may consider the case 7 = const for 
a moment, e.g. 7 = 1 {H = 1/2), corresponding to uncorrelated data. Fig.4 shows the edge 
characteristics of h{q) presented for two distinct series of length L = 2^^, 2^°, generated with 
the autocorrelation exponent 7, and then shuffled. The dependence on 7 is evidently absent 
proving that shuffling procedure was effective enough, while the residual dependence Ah{±15) 
on L is still kept and obvious. 

The detailed analysis of the latter dependence is revealed in Fig. 5 collecting results for 
various data lengths. Astonishingly, this figure suggests a power law dependence between 
Ahi = A/i(7 = 1) and L. 



A/ii(L) = CiL-^S (6) 

where Ci and rji are constant. 

The knowledge of 95% confidence level for this relation is crucial in practise. Its meaning 
is that any result measured above the particular value has probability less than 5%. To obtain 
this confidence level one has to correct Ci and rji parameters by the corresponding quantiles 
calculated from the la uncertainties ac, o"^ of the fit and from the standard deviation 5* resulting 
from the series statistic^: 

Ahf^^iL) = C, exp(/(ac + S))L-^^+f--. (7) 

where / = 1.65 is the respective factor for the particular 95% confidence level. 

Let us take now a closer look at the case of autocorrelated (0 < 7 < 1) finite signals. The 
edge values for h{±15) versus the autocorrelation exponent value 7 were investigated, keeping L 
fixed. Examples of this dependence for L = 2^^ and L = 2^° are shown in Fig. 6. We found that 
cases for other lengths (not shown) look similarly and indicate the excellent linear decreasing 
function of /i(±15) versus 7 in the whole range of autocorrelation exponent. Thus one gets: 

Ah{j,L) = A{L)j + B{L) (8) 

where the coefficients A{L) and B{L) depend on L only. They can be further specified if the 
form of Ahi{L) and Aho{L) = Ah{'~f = 0, L) functions are used as boundary conditions. 

The first boundary condition, i.e. Ahi{L), was already specified in Eq.(6). The profile of 
the second one (Aho^L)) can be deduced from Figs. 6,7. The extrapolation of the fitting lines 
h{±15) versus 7 to the point 7 — )■ (see Fig. 6) gives the collection of Ah{0, L) values, plotted 

^exponential dependence in this formula comes from the uncertainty of regression fit in logarithmic scale 
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against the length of time series in Fig. 7. This can be done for the central values as well as for 

the data satisfying 95% confidence level. It is seen from the Fig. 7 that for fully autocorrelated 
time series (7 0) Aho{L) is represented again by the power law: 

Ah{0, L) = CoL-''° (9) 

with some constants Cq and rj^ to be determined from the fit. 

Linking the shape of boundary conditions given by Eqs. (6) and (9) with the general linear 
dependence in Eq.(8), one arrives with the final formula for A/i(7, L): 

Ah{^, L) = CiL-^'^ + CoL-''°(l - 7) (10) 

The shape of the 95% confidence level for multifractal background noise will be given by 
the same formula but with different coefficients calculated according to formulas like in Eq.(7). 
The final values of these coefficients are collected in Table 1. 
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Co 


^0 




„95% 
'/l 


(-795% 


„95% 

'/o 


0.603 


0.175 


0.453 


0.124 


0.631 


0.171 


0.484 


0.120 



Table 1: The collected results for coefficients of the fit in Eq.(lO) describing the multifractal 
noise thresholds. Ensemble of 10^ independent realizations of time series was considered. 

Our results may also be presented graphically in a form of 'phase-hke' diagrams (see 
Fig.8) .Three separable areas in (A/i, 7) plane can be distinguished for every L. The first area 
corresponds to multifractality connected entirely with finite size effects. It is marked in red in 
Fig.8. The second domain, marked in light green, is related to {Ah, 7) range where multifrac- 
tality may occur due to the long memory present in data but independent on the chosen time 
scale. The 'true' multifractality, i.e. the one related with long memory entirely dependent on 
the time scale, may occur only in the white region (at 95% confidence level). 

3 Multifractal spectrum analysis of finite size effects. 

The multifractal spectrum width is considered as another useful measure of multifractality 
included in analyzed signal. Analogically to the Ah{L, 7) analysis, one may ask for the depen- 
dence of multifractal spectrum width Aa on the signal length L and on its persistency level 
7. These results are obtained with use of Eq.(2) applied to previously discussed generalized 
Hurst exponent calculations. Thus, the results should lead to similar qualitative conclusions, 
nevertheless quantitatively they might be also valuable from practical point of view. 

The examples of multifractal spectrum [a, /(ct)) for finite monofractal signals are shown in 
Fig.9. Three cases: for strongly autocorrelated, medium autocorrelated, and weakly autocor- 
related signals are considered there for two distinct lengths of data: L — 2^^ and L — 2^°. 
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As previously, the first step is to examine the Aa characteristics obtained for randomly 
shuffled signal (7 = 1). Fig. 10 shows the minimal a^m and maximal oimax value of a param- 
eter revealing lack of dependence on 7. This proves again that shuffling was sufflcient. The 
dependence of Aa(7 = 1) = Aai on L obeys a power-law relation: 

Aai(L) = (11) 

shown in Fig. 11, with unknown constants and i^i. The 95% confidence level is given by 
equation analogical to Eq.(7). 

The edge values for aminimax as a function of 7 are presented in Fig. 12 for particular lengths 
L = 2^^ and L = 2^*^, to confront them with previous findings for /i(±15). The linear dependence 
for all other lengths (not presented) is also observed. Once we repeat the same approach as in 
the previous section to Aa(L, 7) dependence taking into account the second boundary condition 
AaQ{L) = Aa{L, 7 = 0), we arrive with the final formula describing the character of multifractal 
spectrum width (see Fig. 13), similar to the one found in Eq.(lO): 

Aa(7,L) =DiL-«i7 + Do^-«°(l-7)- (12) 
The values of fitted parameters are gathered in Table 2. 





ei 


Do 


eo 


^95% 


^95% 


n95% 


f95% 
SO 


0.686 


0.129 


0.572 


0.089 


0.784 


0.120 


0.670 


0.079 



Table 2: The results of the fit for coefflcients in Eq.(12) done on ensemble of 10^ independent 
realizations. 



4 Concluding remarks 

We have shown qualitatively and quantitatively how multifractality arises from the finite size 
effects and (or) from autocorrelations not changing with the time scale being formed by the 
specific autocorrelation exponent 7. This kind of multifractality, called by us 'multifractal 
noise', should be clearly distinguished from the 'real multifractality' caused by memory effects 
dependent on the time scale and thus leading to different scaling properties at various scales. 
We provided analytical formulas describing the multifractal noise threshold which turns out to 
be the power law function of time series length L and the linear function of autocorrelation ex- 
ponent 7. Our approach differs from the one presented in Ref . |19] where long memory in signals 
was produced from the induced power spectrum profile instead of the direct autocorrelation 
input between data resulting immediately from Eq.(3). 

Two description methods for multifractality were considered, i.e. the generalized Hurst 
exponents and the multifractal spectrum analysis. In both cases multifractal residual effect in 
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monofractal finite signals has been found. We have shown that this effect measured by the 
spread of generahzed Hurst exponent or the width of multifractal spectrum grows hnearly with 
autocorrelations level in time series and decays according to power-law with their length. We 
have estimated numerically the level of such multifractal noise and we captured it in simple 
analytical equations. 

Finally, one should compare the obtained multifractal noise threshold with examples of 
the real multifractal data. We took them from finance because of common agreement that 
multifractality is a characteristic feature of financial markets. This problem is considered in 
Fig. 14, where the simulated 'phase-like' diagram for data length L ~ 2 x 10^ is shown, together 
with multifractal properties of various markets - both for price indices and for volatilities. 
The particular length L has been chosen as the average length of available data for analyzed 
markets. The multifractal properties for price indices were taken from ^7] , once the respective 
features for volatilities were originally calculated by us for the purpose of this paper and are 
based on historical data available in web [28j. 

The multifractality for the volatility series is noticed. However, the presence of multifrac- 
tality for the price index data is already not so obvious for all markets. One may find indices 
where multifractality comes indeed as a result of scaling properties changing with the time scale 
(e.g. Venezuela, Indonesia, China). Simultaneously, there are markets where the observed mul- 
tifractality is generated mainly (e.g. Philippines, Taiwan, Germany) or even entirely (Ireland) 
by the finite size effects. In the case of Philippines, Taiwan, Thailand, Germany, Spain and 
Greece almost 80% of multifractality in price indices is related to such an effect. 

Thus, the multifractal properties of real financial data, may substantially (even in 80% — 
100%) originate from the multifractal noise, what makes difficult in some cases to separate 
what main phenomenon is really responsible for the effects one observes. This confirms that 
multifractality is very tiny and delicate effect and one should be especially careful drawing far- 
reaching conclusions from the multifractal analysis in finance and in other areas. Our formulas 
are general enough to be applied also to other kind of real data in order to distinguish if and 
how their multifractal properties result from the pure multifractal origin. 
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Figure 1: Efficiency of FFM for replication of autocorrelation properties in time series. The 
examples for input values 7 = 0.2, 0.5, 0.8 are shown in log-log scale for the generated data 
of length L = 2^°. The lines present the fit to the desired power law dependence of Eq.(3), 
while error-bars show la standard deviation following from the considered statistics of 10^ 
independent realizations. The output 7 values from the fit are found 'jout = 0.203(±0.009), 
0.498(±0.012), 0.782(±0.054) respectively. 
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Figure 2: Scaling of g-deformed fluctuations within MF-DFA. Results are presented for two 
different lengths of time series L = 2^^, 2^*^, three autocorrelation parameters 7 = 0.1, 0.5, 
0.9 and q = —15, —10, —5, 0, +5, +10, +15. All plots confirm the proposed scaling range from 
r = 10 till r = L/4. 
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Figure 3: Generalized Hurst exponent for monofractal signals generated within FFM with 
various autocorrelation properties. Two cases, for L = 2^^ and L = 2^*^ are shown with different 
autocorrelation level. Error bars correspond to statistics of 10^ generated series. 
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Figure 4: Edge values of the generalized Hurst exponents h{q) for two different lengths of time 
series L = 2^^, 2^" constructed with long memory present (7 < 1) and then shuffled to kill this 
memory. Dependence on the data length is readable. 
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Figure 5: Spread Ah of generalized Hurst exponent versus length of data drawn in logarithmic 
scale for the signal with no memory. Power-law dependence between Ah and the data length 
is visible. Data point corresponding to L = 2^ (marked as triangle) is slightly above the fitting 
line due to insufficient statistics for so short signal. Therefore, this point has been neglected 
during fitting procedure. Results of the fit are Ci = 0.603, rji = 0.175 for the central values 
and Cf^'^" = 0.631, r]f^'^'' = 0.171 for 95% confidence level. The latter fit is marked as the blue 
top line. 
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Figure 6: Edge values of generalized Hurst exponent for series with long-memory. Figures show 
the linear dependence between the edge values /i^ = /i(±15) and 7 exponent. Extrapolation of 
fitted lines to the point 7 = is interpreted as the edge values for fully autocorrelated signal 
{C{r) ^ 1, Vr) 
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Figure 7: Spread Aho of generalized Hurst exponent for fully autocorrelated time series (7 = 0) 
versus the length of data. Green line presents the power-law fit in log scale and the blue line 
corresponds to 95% confidence level resulting from statistics. Fitted parameters are Cq = 0.453, 
T]o = 0.124 and C^^^° = 0.484, r^^^^" = 0.120. Data point corresponding to L = 2^ has been 
removed from the fit due to insufficient statistics for so short signal leading to huge uncertainty 
in the estimation of generalized Hurst exponents within MF-DFA. 
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Figure 8: Examples of phase-like diagrams for two data lengths L = 2^^ and L = 2^°. Three 
separable areas can be distinguished, marked in red, light green and white. They correspond to 
{Ah, 7) domains where multifractahty is caused at 95% confidence level respectively by finite 
size effects, constant autocorrelation level independent on the time scale and 'true' multifrac- 
tahty (white). 
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Figure 9: Examples of multifractal spectra for finite monofractal strongly autocorrelated (7 = 
0.1), medium autocorrelated (7 = 0.5) and weakly autocorrelated (7 = 0.9) signals. Mean 
results for two distinct lengths of signal L = '}P and L = 2^° are compared for statistics of 10^ 



series. 
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Figure 10: Edge values of the Holder parameter ot for two different lengths of time series 
L = 2^^, 2^° constructed with long memory present (7 < 1) and then shuffled to kill this 
memory. Dependence on the data length is readable. 
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Figure 11: Spread Aai of multifractal spectrum versus length of data drawn in logarithmic 
scale for the signal with no memory. Power-law dependence between Aai and the data length is 
confirmed. Results of the fit are Di = 0.686, ^1 = 0.129 for the central values and ^° = 0.784, 
^95% ^ Q^i20 for 95% confidence level. The same notation applies as in Fig.5. 
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Figure 12: Edge values of Holder parameter a for series with long-memory. Figures clearly 
show the linear dependence between the edge values aminimax and 7 exponent. 
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Figure 13: Spread Aao of multifractal spectrum for fully autocorrelated time series (7 = 0) 
versus the length of data. Fitted parameters are Dq = 0.572, ,^0 = 0.089 and Dq^^" = 0.670, 
^95% _ 0.079. Notation is inherited from Fig. 7. 
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Figure 14: Simulated 'phase-like' diagram for L = 2 x 10^ with market data charted. Green 
squares represent multifractal properties of price indices and blue dots represent multifractal 
properties for volatilities. The Bloomberg code has been used to describe markets. The source 
of historical data is indicated in text. The corresponding values for Hurst exponents H are also 
indicated (top axis) for convenience. 
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